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Abstract 



We suggest a new asymptotic representation for the solutions to the 2-D wave e 
quation with variable velocity with localized initial data. This representation is a gen- 
eralization of the Maslov canonical operator and gives the formulas for the relationship 
between initial localized perturbations and wave profiles near the wave fronts includ- 
ing the neighborhood of backtracking (focal or turning) and selfintersection points. We 
apply these formulas to the problem of a propagation of tsunami waves in the frame of 
so-called piston model. Finally we suggest the fast asymptotically-numerical algorithm 
for simulation of tsunami wave over nonuniform bottom. In this first part we present 
the final formulas and some geometrical construction. The proofs concerning analytical 
calculations will be done in the second part. 



1 Main equations and a simple example: the 
wave field in the case of constant bottom 

X 

1.1 Some notation 

Let us introduce the notations used in this paper. A two dimensional vector can be 
written with capital or small letters X = (Xi,X2) or x = (x\,X2). The vector can 

be written also as a column vector f ^r^) ' Two dimensional vectors X and Y can 

form a column vector ^y^j with four rows. The real scalar product between two 

vectors X and Y, with real components, is indicated by < X, Y >, the complex scalar 
product among bi-dimensional vectors Z, W, with complex components, is written as 
< Z, W > c , the two by two matrix generated by two bi-dimensional vectors X, Y is 
written as (X, Y) where in the first column there are the components of the vector X 
and in the second column those of the vector Y; the transposed matrix of C is denoted 
by t C. 
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1.2 Main equations 

Let us remind the statements of problems used in tsunami wave problems (PQ-^O]) as 
well as in general linear water wave theory (|ll|-[T7j. 

Let us assume that the bottom of the basin is moving H = Hq(x) + Hi(x,t). 
We assume also that the perturbation H\(x,t) is small with respect to Hq 
\H\\ « Hq(x), and that Hi is localized in a neighborhood of some given point xq. If 
L is the dimension of the region where the wave phenomena is studied, and I is the 
dimension of the perturbed region, then our hypothesis implies that I « L. Another 
assumption is that the bottom "changes slowly", i.e. that VHq ~ /i, where \x is some 
small ("adiabatic parameter"). We discuss below its meaning. Introducing the scaled 
variables x' = £ , then H = H (x') + xH\(^,t), where fj, = j << 1. 

The equation for the velocity potential $ in the water —H < z <rj, where rj(x, t) is 
the sea elevation in the linear approximation, has the form, in dimensional variables,: 

A$ = 0, (1.1) 

Vt- lr z=o = 0, (1.2) 
Oz 

*t + OT|*=o = 0, (1.3) 

9$ <9<£ 

— = — + < VH, V$ >= v(x, t) z =-H- 1.4) 
on oz 

where v(x,t) is the normal component of the velocity of the motion of the bot- 
tom in the point x. The velocity v can be expressed by means of the derivative 
by : / \J (VH) 2 + 1, since V is the projection of the velocity on the vector 

1 =*(Vij", 1) normal to the surface z = —H. If we consider X7Hn to be small (be- 

V(v/?) 2 +i v 5 ; u v 

cause of the slow variation of the bottom relief), and that also Vi?i is small (because 
of the small amplitude H\), then we have 

dHi 

1.3 The solution in the form of the Fourier transform 

Let us begin considering the system ([1.1)1 - (|1.4|) in the case of constant bottom. In this 
case the velocity potential and its derivatives are zero for t = 0. We make the Fourier 
transform of the system ()l.l|) - (|1.4j) with respect to the variables X\,X2- The dual 
variables will be denoted with pi,P2 and the Fourier transform of the corresponding 
function will be considered as a "wave". Then (|l.lj) - (|1.4)l gets the form 

$^-p 2 $ = 0, (1.5) 



»7t - -a-U=o = 0, (1.6) 
oz 



($ t + g?f)\ z=o = 0, (1.7) 



d<S> . dH-L 
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Solving lE^D-JTEJ), we find 



and 



~_ cosh(( Z + g)|p|) sh(z\p\) djh 

cos\iH\p\ \p\ cosh.(H\p\) dt 



~ 1 dH 



Thus the equations (|2.8[1 - (|3.2|) take the form 

^ - \p\ tanh(H\p\)ip } tt . , = 

dt W] y inj ^ cosh(iJ|p|) dt 

^ + OT = (1.11) 

Where (p = & t \ z= o, and we have the initial conditions t = 

(p\t= Q = 0, (p t \t=o = 77U=o = 0. (1.12) 



These conditions define the so called Cauchy-Poisson problem for the system 1)1.11(1 . 
They are compatible with the perturbation of the bottom only if we suppose that the 
earthquake starts at at a time different from zero. So we assume that the bottom has 
an "instantaneous" movement at a small time t = e: 

H x (x,t) = 6(t-£)V(x), (1.13) 

then we send e to zero at the end of the calculation; the smooth function V(x) decays 
rapidly at infinity. 

Differentiating the first equation in 1)1. 11|) with respect to t and substituting §j| 
with —gfj we get the equation for fj: 

Vu + Crj - = 0, C = 9 \p\ tanh(tf M). (1.14) 

Differentiating the second equation of the system ()l.llj) with respect to t and sub- 
stituting the derivative rjt with the expression of the first equation and consider- 
ing the condition that the source is active at the moment t = e > 0, we get 
¥?tt|t=o = —g\p\ tanh(H\p\)ip\ t= o = and the initial condition for 1)1.14)) 

Vt=o = Vt\t=o = 0- (1-15) 



It is easy to find the solution G of the homogeneous equation associated with ()1.14)) : 
G u + C(p, H)G = 0, G\ t=T = 0, G t \t=T = 1, 

G{t ^ p) ~ 2VZ " 7z • 

In this way the solution of the non homogeneous equation 1)1.14)1 is 

- fnu 1 d 2 H x {T,p) j 

V = J G{t ^ P) cosHH\p\) dt- ^ 
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The inverse Fourier transform of the function fj gives the elevation of the free surface. 

Hi(t 



Under our assumption of instantaneous motion at time e we have 9 H f)}i' P ^ = S'(t — e)V 



and so: 

V d f s3nV£(t-T) s , V 



coshH\p\dr ^[C_ ' lr — e cosh^T|p| 



cos 



We send now s to zero so we get the function fj = cosh ^j-| p |) cos y/C(t). It is evident 
that fj is the solution of the equation Q1.14|) with the following Cauchy conditions 

f}\t=o = — rjw\ iv fj'\ t =a = 0. (1.16) 
cosh(ii |p|J 

We shall discuss the meaning of such initial conditions for the function r/ later. 

1.4 Solution of the Cauchy problem for constant bottom 
and instantaneous source 

Let us study the solution r\ corresponding to (|1.16j) . It is not restrictive to assume that 
the center of the source is located in the origin of the coordinates xq = and that 
the perturbation decays rapidly with the distance from the origin and that it has a 
maximum in a small neighborhood of the origin. We use also dimensionless variables: 

where I is the size of the shifted region and 

V = h I V&e-^dt = ±f V{y)e- a <™>dy = IV (pi), 

i i 

cosh(\p\H) 



where we made the substitution (y = 4, £ = yf) and V(p) is the usual Fourier 
transform of the function V(y). We assume that V(y) is smooth function rapidly 
decaying as \y\ —> oo. 

Then we can make the inverse Fourier transform: 



V f e ±u ^ c ^ +i<p ' x> fi (p)dp = — S± / e ±u V c (P' H )+^> }—^ V (pl)dp. 

J 4ir J cosh(\p\H) 

Changing the variables pi = p',p = p'/l, we get 



V 4vr , 



V = / e V i ufi , , — V(p)dp. 

47T J cosh(|p|-j-J 

In this way the problem is reduced to the computation of the asymptotic behavior of 
the integral. 

We will study the asymptotic values for \x\ » I. We change to polar coordinates 
(p,ip) in the integral, where 93 is defined as the angle among p and x — xq. Thus 
p = p0(</?)#r, where &(<p) is the two dimensional matrix defining the rotation of an 

\x\ 

angle ip. 

„ , , / cos ip — sin 03 N 
9(v3) = . 

Vsmy? cos ip 
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Then the last integral has the form 



4vr 



2tt 



pdp I dip exp 



± ityj^j- tanh(p— )) exp 



cosh(pe^f) 



The internal integral can be computed using the method of stationary phase. The phase 
has the form: $ = cos p, the equation = gives p> = 0, (p = it; furthermore 
it is not possible to apply the method of the stationary phase in the point p = 0. 
Nevertheless one can take a sufficiently small neighborhood of the saddle points of the 
variable p and show that, 1191 12UI I24j . it is smaller than the contribution of the 
terms that we neglect. The result is: 



if? I a: I 

p / dpe^r 1 ™^ 
Jo 



2tt 



cosh (P^f) 



cosh(pe^f) |x| 




I / '<> ie\A~ X . j /4 / lg - , 



and 




■'I 7o cosh(p T :c " 



R i 



exp ( ± ity — tanh(p— (e 



e i 



n±Pn))- 



Let us consider the last integral. Its complete phases are: 



<T>±.±/7 - - ; (jl p\;\n\\{j)—j-) I p x ; 7. 



For i > 0, p > the derivative 



<9p 



is strictly positive, this implies the absence 



of critical points for the functions 3>±,+. It follows that these terms give, for t > 0, a 
contribution to the wave field which is asymptotically small with respect to the other 
contributions and so it can be dropped. Furthermore since V is a real function then 
V(p-nj) and V{— Pjfr) are complex conjugates so the last integral may be written in 
the form 




Re / dp- 

lo cosh(p|f rT 



X/ Wv(^)e-W4 exp (| (eW _ t/3 ^H J-L tanh(^)) . 



Since the ratio y- is rather small, the source is localized, and the function 
decays rapidly as a function of p, then the main contribution to the last integral is 



coming from the small values of p. Then we get that the functions 



i 



cosh (p j 



T) 



and 



tp^gHy -Jjj tanh( £ p) can be expanded in Taylor series. If we substitute the first 
function with 1 we neglect a term of the order of O(-j-) 2 - The second function can 
be approximated by the first two non zero terms of its expansion tp^/gH(j — \{^j-) 2 ) 
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making an error of the order of t^/gH^) 4 - It is clear from the previous estimates that 
these terms are small and so we obtain 

oo ■ 2 

jMp-W4_ n ^ n |„| _tm/ffHYl - f —(-^ 



d^^V(P n )e- W4 exp (-( e |x| - WgH(l - ^( T ) 2 ))). 



X | y e v r cl 1 r v ° v 6 1 

It will be explained below that the integral gets its larger values in the neighborhood 
of the front, i.e. near the curve ( circle) \x\ = y/gHt. In this way the dispersion effects 
can have an influence on the asymptotic values in the far wave field under the condition 
that the coefficient of p 3 in the exponent is larger or equal to one. Thus we obtain 
different behaviors, putting \JgHt equal to |x| in this coefficient, according to different 
relations among \x\,H,l (compare [I]-|lUj.[21]): 

a) For |ar| >> -jp the dispersion has an important influence in the neighborhood 
of the front, and the asymptotic can be expressed by means of a function similar to 
the Airy function. In this case the behavior of the function V is not important for the 
definition of the profile of the front. 

b) For \x\ ~ jp the weak dispersion and the function V have equal influence on 
the formation of the wave profile; 

c) For | 

x \ << jj7 the dispersion is not important. If the term with p 3 , is dropped 
from the phase of the integral an error of the order of \x\H 2 /I 3 is done. 

Let us consider the example where H = 4km, I = 40km, thus I 3 /H 2 = 4000/cm. 
Thus a (weak) effect of the dispersion starts at 4000/cm. If the size of the source 
increases twice this distance increases 8 times and becomes 32000 km, a distance larger 
than any ocean. Thus we will start analyzing the point c (it possible to neglect the 
effect of the dispersion). 

1.5 Asymptotic behavior of the wave field with very small 
dispersion in the case of constant depth 

NEW 

Thus, assuming that the inequality \x\ « I 3 /H 2 is satisfied, we have 

r°° ~ x / i , 

Re dp^Y(p—)e- m / 4 exp - tpy/gt) = 

4^Re(e-/ 4 F(^,ii)), $ (x , t ) = |x| -t^gH, 



1 I 



V2ttV \ x \ 



where 



If 00 - r 

F(z,n) = — / e^^pV(pn)dp n = — . (1.17) 

V Z7T Jo \X\ 

Here n is a unit vector. It is natural to introduce its angle ip in such a way that tp = 
corresponds to the axis x\. Then 

"="»> = ("$)■ <" 8 > 

Hence the function V(p, n(ip)) depends on (p,ip) and the function F(z,n(ip)) depends 
on (z, tp). Not to complicate notation we use the same symbols V and F for them and 
sometimes write V(p,ifi) and F(z,ip) instead V(p,n(ip)) and F(z, n(ip)) respectively. 
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We note, that the function F(z,n) decreases for \z\ — ► oo as an inverse power. 

2 

Indeed, let us change variable in the last integral p = \\ then 

1 f°° 2 ~ V 2 

n) = -=Re / y 2 e iz < V"f ) V(^n)dy. 

Let us use the method of the stationary phase, we get, because of the presence 
of the factor y 2 under the integral, F(z,uj) ~ > ^ ^(0) 7^ 0- Thus for 
\\x\ — y/gHt\ >> I and \x\ >> I, we have that n ~ nj7(0). 

II \x\ 

Example 1. Let us give some example of the function F(z,to). We choose for the 
function V , defining the source, the function 

V{y) = V cos (aiYi + a 2 Y 2 + x ) e - h ^-^ , Y = 0(%, (1.19) 

cos 6 sin ( 



9(0) 

sin t> cos ( 

where V ,ai,a 2 ,b\, b 2 > 0, 9,x o-re parameters. In this case the function F(z,ip) can 
be expressed in terms of parabolic cylinder functions D_ 3 / 2 or confluent hypergeometric 
functions iF x (see J2l\/ . 3462, page. 351) 

V(p, = ^§=e- a -^ cosh(*<5 + jp), (1.20) 



FM ) = *l Re ( e -f r v ^(e-4 a -HwH^e w + e-4 a -Tr^» e - w )dp) = 
2y2n Jo 

VVhh ( ^ (7 + iz) 2 7 + i0 (-7 + iz) 2 + - W \ _ 

= 4^7l Re(exp( ^^ )D - 3 /^^ )e +exp( ^^ )D - 3/2( "^~ )e ) = 

^vS fe[(Q++Q - )] ' 

^(^) = ^^(^r(i T i) ^4 I, f) + ^r ( i) ^4 §, d)), 

where T is a gamma function, and 

w± = ±j + iz, a = (biah 1 + b 2 alj) / (4b\b 2 ) , 
(5 = (6isin 2 (^ - d) + 6 2 cos 2 (V - 0))/(4M 2 ), 
7 = (&ia2 sin('0 - 6) + 6 2 ai cos(V> - 0))/ (2b\b 2 ), 

The figures Fig. 1, Fig. 2 present some types of profiles. 

Main conclusion: the phase in the neighborhood of the front defines completely 
a one parameter family of trajectories which generate the front. Further we remark 
that, since the function F decreases in a neighborhood of the front, we can expand in 
the formula IJ1.17J) |x| in a neighborhood of the front, keeping in the expansion only 
the zero order term, and that we can substitute the factor — }= ( the amplitude of the 



wave) with the term , We want to find analogous formulae for the wave field in 
the case of negligible small dispersion and for variable bottom. 
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2 Asymptotic behavior of the wave field over 
nonuniform bottom for very small dispersion 

2.1 The wave equation, rays and wave fronts 

In this section we start the analysis of the behavior of the amplitude of the wave when 
the bottom is not constant. We use here well known objects and their characteris- 
tics which one can find in books connected with the semiclassical asymptotics and 
ray method, geometrical optics and wave fronts, Hamiltonian mechanics, catastrophe 
theory etc. We try to collect here all necessary objects and give their description in 
elementary form. More complete form and details one can find in |25|-|51| It is clear 
that in practice we have studied the solution of the wave equation in the previous sec- 
tion. In order to be accurate we introduce in the calculations the characteristic depth 
Hq of the basin, the small parameter 

M = { (2-1) 

expressing the relationship among the characteristic size of the source and the charac- 
teristic size of the basin. We begin introducing non dimensional variables in the equa- 
tions. Then after a suitable change of variables x = x' /L, t = t 1 \J gHo/L, H = HqH'(x') 
our equations and initial data will take the form: 

^=g<V,H(x)V> V , (2.2) 
V\t=0 = V{-), r}t\ t=0 = 0. (2.3) 

Our asymptotic expansions will be done in term of this parameter under assumption 
[i -C 1 and that the domain in and time interval [0, t] G while asymptotic 
expansion are working independent of p. To come back to original variables it is enough 
to change in final asymptotic formulas fj, by I. 

We assume that the source of the perturbation is localized in x = 0. It is easy to see 
that finding the field far from the source,|x| ^> I, is similar to find the asymptotic values 
for (i — > in the problem (|2.2|) . The problem now is to study the wave equation with 
variable coefficient. The asymptotic values of the wave amplitude r] can be expressed by 
means of the wave front formed by rays. It is a known fact that instead of the straight 
rays one has to introduce curved rays and characteristics which are the one dimensional 
family of trajectories P(ip, t),X(ip, t) of an appropriate Hamiltonian system . The ends 
of the rays again form the wavefront, but now it can be different from the circle, a 
more complicated closed curve probably with cusps and self intersection points. In the 
considered situation these rays and characteristics are determined in the following way. 

We introduce the function C(x) = gH(x) and, as before, let n be the unit vector 
Q1.18JI directed as the external normal to the unit circle. Then the Hamilton system is: 

p dC 
x = -^C(x), p = -\p\ — , x\ t=0 = 0, p\ t=0 = n(ip), (2.4) 

i.e. the family of trajectories P(t, ^),X(t, ip) going out from the point x = with unit 
impulse p = n(ip). Let us indicate C(0) = Cq. The Hamiltonian corresponding to ()2.4|) 
is 7i = C{X)\p\. From the conservation of the Hamiltonian on the trajectories we have 
the important equation 

\P\C(X) = C . (2.5) 



S 



The projections x = X(t,tp) of trajectories on the plane are called the rays. Recall 
that the front in the plane M 2 at the time t > is the curve 74 = {x 6 K 2 |x = X(ip, t)}, 
The points on this curve are parameterized by the angle ifi £ (0, 2ir]. If in each point 
x of the front 7^ ^ 7^ 0, then the front is a smooth curve. The points where 

S = are named /oca/ (backtracking or turning points), in these points the front 
looses its smoothness. In the situation in which the focal points appear, (they are very 
interesting from the point of view of tsunami), it is reasonable to introduce the concept 
of the front in the phase space x at the moment t > 0, i.e. the curve Tt = {p = 
P(ip,t),x = X (tp , t) , ifi 6 [0, 27r]}. We note that at least one of the component of the 
vector P^,X-q, is different from zero, see Lemma 3; and also the rays x = X(t,ip) are 
orthogonal to the front jt'- (X, X^) = see Lemma 3. 



2.2 The wave field before critical times. 

It is not difficult to check that a (possibly sufficiently small) t\ exists such that, for any 
t, t\ > t > 5 > 0, there are no focal points in 74. The first instant of time t cr , in which 
focal points are formed is called critical. Let us first write the solution before critical 
times, larger than 6, when the front is already defined. In this case the asymptotic 
solution is defined in the following way. We define a neighborhood of the front for 
sufficiently small ( but independent of h) coordinates ip,y, where \y\ is the distance 
among the point x belonging to a neighborhood of the front and the front. For this aim 
we will take y > for the external subset of the front and y < and for the internal 
subset of the front. Then a point x of the neighborhood of the front is characterized 
by two coordinates: ip(t,x) and y(t,x), where tp(t,x) is defined by the condition that 
the vector y = x — X(ip, t) is orthogonal to the vector tangent to the front in the point 
X(ip,t). Thus we have the condition {y,X^{ij),t)) = 0. Let us find the phase 



S (t, x) = (pm, «), x - x m , x), t » = C(x ™ xU)) v = ^ H{X "Z), t )f 

The second equality is a consequence of the equation (|2.5j) . 

Now we state the first important theorem of this paper connecting the wave ampli- 
tude with the initial perturbation V(x) and the profile of the bottom and the integration 
over the characteristics. 

Theorem 1. For t cr > t > 5 > in some neighborhood of the front jf, not depending 
on \i, T], the asymptotic elevation of the free surface, has the form: 



V " " m Re[ e -fF(^),„(rt)]L_, (til+ 0(^). (2.6) 



Outside this region r\ = 0(fi 3/2 ). The function F(z, n) is defined in p.!7jl . 

In this way till the critical time the asymptotic elevation of the free surface is 
completely defined by means of the trajectory, which forms the front of the wave, and 
of the function V, corresponding to the source of the perturbation. Despite of the 
simple and natural form of the asymptotic of r/, the proof of the formula ()2.6|) is not 
trivial at all; the main step is the computation of the function V, more exactly the 
proof of the fact that the formula is the same as in the case of constant bottom, if the 
right choice of the rays is made. We will give below a constructive approach of the 
proof of this formula, in the meantime we now show some elementary consequence of 
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the equation (|2.6I) . Since the phase S(x,t) is equal to zero on the front and S(x,t)//i 
gets large going out from the front, then 77, as one could expect, decreases enough 
quickly and the maximum of \rj\ is attained in a neighborhood of the front. As a 
consequence, rj can have some oscillations depending on the form of the source. The 
second factor in (|2.6j) is the two dimensional analogue of the Green law, well known in 
the theory of water waves in the channels: the amplitude rj increases when the depth 
decreases as the inverse of the fourth root of the depth 1 / ^C{x) = 1/ \j H(x); the 
factor l/y^A^I is connected to the divergence of the rays, in other words if a smaller 
number of rays goes through a neighborhood of the point X(ip,t), the smaller will be 
the amplitude of the wave field. The factor C ( X (fp°t x) t)) appearing in the formula of 
the phase expresses the phenomena, also well known, of the "contraction" of the wave 
profile as the depth decreases and the increase of its amplitude. In fact the amplitude 
increases because of the factor in front of the function V but also the phase S(x,t) 
increases and this makes the wave profile narrower. This result explains the well know 
fact that the wave length of the tsunami decreases when the wave approaches the coast 
and that its amplitude increases. The same profile (i.e. a section of rj(x, t) for fixed t 
and ip) can depend on the way the trajectory ( ray) intersects the initial perturbation 
of the bottom at t = 0. It is just this fact to give the dependence of the diagram of 
the directions on two factors: the shape of the source and the angle of its intersection 
with the ray passing through a given point of the front. For this reason, depending on 
the form of the bottom, two rays going out with two very different angles, can arrive 
near a point of the front and contribute to the profile with very different amplitudes. 
These effects can be well seen in the figure Fig. 3 

Remark 1. The main argument usually used for deriving the the analytical and and 
asymptotic formula of the solution of (j2.2j) is that, after the front is formed, the prob- 
lem is essentially one dimensional in space till the appearance of focal points and its 
dynamic is described by the one dimensional wave equation with non uniform velocity. 
But it is possible to construct different particular wave solution, localized in different 
neighborhoods of the same front but having completely different profiles. The question 
is the right choice of the function describing this profile, i.e. a function connected 
already with the construction of the solution of bi- dimensional problems, containing 
mathematical difficulties such as the presence in the solutions of the wave equation 
( with variable coefficients) of effects of intersection of the characteristics (which in 
quantum mechanics are also called terms). This happens for very small wave vectors, 
i.e. only for very long wave lengths. A general approach allowing to treat difficulties of 
this type were developed in \21j/ and, in particular for the Cauchy problem with localized 
initial conditions for the wave equation and hyperbolic systems, in the papers \20\ 
These approach is basically founded on the following argument. The solution has two 
contributions, the first, corresponding to "very long" waves, can be found only directly 
by numerical methods. It has a very small amplitude and, in the case of the problem of 
the tsunami does not have a very great importance. The second has a wave length long 
compared with the depth of the basin of the wave but sufficiently small compared with 
oceanic scale, in such a way that it is possible to apply effective asymptotic methods 
as the expansion in rays. These arguments for the given problem ()2.2|) were given ac- 
curately in the papers \22\ \2ty . but the final formulas, based on the asymptotic for the 
system with constant coefficients \47j , are not very efficient from the point of view of 
practical applications. The basic of the derivation of the formula ()2.6|) of [411 consists 
of two steps: 1 ) the construction of the asymptotic expansion in the smoothness of the 
fundamental solutions of the problem (J2.2I) ( parametrix) and 2) the evaluation of the 
asymptotic of the convolution of this solution with the initial function 770 . Basically 
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our main observation ( missing in the papers \4 r Fj ) is that this last asymptotic may 
represented in a simple and useful way (|4.1|) . from which the representation Q2.6JI . as 
well as representations (|2.9|) . (|2.1(ij) working after appearance of critical points follow 
immediately. 

2.3 The structure and metamorphosis of wave profiles. 

2.3.1 The Maslov index and metamorphosis of the wave profile. 

For t > t CT when the focal points appear, as it is well known in the wave theory, the 
front can have the "angles" and sometimes the front lines can have self intersection 
points. The ends of the arcs corresponding to these angles are the focal points (or 
backtracking or turning points). For t > t CT the front divides in some arcs 7^, indexed 
by the number j, separated by focal points. The internal points of these arcs are the 
ends of the trajectories P(ip, t), X(ip, t) with the same topological structure. Namely 
these equivalent trajectories cross the same numbers of focal points at times t F before 
t, t < t. They are characterized, from the topological point of view, by the Maslov 
index, an integer number m(ip, t) depending on if), t. The Maslov index m can be 
defined on the regular points of the front in different ways, we give below a more 
practical definition of this important concept by a simple definition of its increments 
in the problem under examination. The index m is related to the sign of the Jacobian 
J = = (X,Xy,). The function J is equal to zero in the focal points and only 

in these points. Thus moving along the front jt or along the trajectory (P, X) after 
crossing the focal point, the Jacobian can change its sign. Actually the Maslov index 
prescribes a receipt for assigning the correct sign to the square root of J and it can be 
defined in a way independent from the trajectories. But if we move along a trajectory 
there is, in this problem, the nice and useful fact that the Maslov index coincides with 
the simpler Morse index. So, considering the trajectories arriving to 7^, we have that 
the Morse index m(ijj,t) of the point x = X(ip,t) S M 2 is equal to the number of focal 
points on the trajectory p = P(i(),t), X(i(),t),t G (0, t) arriving to x = X{ip,t). Note 
also that, as the time t changes, the ends of the arcs 7;? produce the entire set of focal 
points. It is also a well known fact that these sets constitute the (space-time) caustics 
which are the singularities of the projections of some Lagrangian manifold (we denote 
it M 2 ) from the phase space 1^ to the plane (configuration space) M 2 . 

Example 2. 

Let us illustrate the concepts explained above by the example [] about the waves on 
an axially symmetrical bank described by the the depth function 

H = H(p),p=^x 2 1 + xl (2.7) 

In this case there exists an additional integral 

P v = X1P2 - x 2 pi (2.8) 

and the Hamiltonian system \2. 7| ) is completely integrable. 

We assume that the source is located in a neighborhood of the point x\ = 0, x<± = 
—po- After the appearance of the focal points one has the following picture (see fig. 4 )For 
each fixed time t the front "ft is separated into two arcs: the first, a long one, is 7* with 
self-intersection, and the second, a short one, is 7^, located between the angles on the 
fronts. The union of the ends of the arc 7 2 for different times t gives a caustic. The 
arc 7 t x consists of the ends of trajectories (rays) without focal points on them (except 
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t = 0). Thus the Jacobian J(ip,t) = detfX, X^)(?/>, r) > for fixed ip and for each 
t £ (0,i]; hence the Morse index m(x £ 7^) = 0. On the contrary the arc y% consists 
of the final points of the trajectories (rays) which cross one focal point at some time 
t = tir(V'),0 < tpfy) < t when they touch the caustic. In this case before tp(ip) J > 0, 
J(tp,t F (tp)) = 0, and J < fort > tp{ip). Hence m(x G t|) = 1. 

Now let us fix the time t and move along the front ■ Then after the passage through 
the focal points the phase — 7r/4 in formula (|2.6f) increases by a quantity — 7r/4 ± 7r/2, 
where ±1 is the jump of the Maslov index. Finally after passing through several focal 
points instead of the factor e 4 one has the factor e 4 2 The number m is 
defined mod 4. The appearance of this new factor produces crucial changes of the form 
of the wave profile in the formula l|2.9|) i.e. in the function Re(e~~ 2~F). This fact 
is analogous to the well known metamorphosis of the discontinuity in the theory of 
hyperbolic systems, and the formula (|2.9|) describes explicitly the appearance of the 
same fact in the case of localized initial perturbations. 

Example 3. Let us give the examples of the transformation (metamorphosis) of the 
wave profile depending on the index m and on the source of gaussian type. We have the 
following pictures for m = 0, 1, 2, 3, mod 4 (see fig ????). Thus in the cases m = 2, 3 
we have the "overturned" profiles with respect to the cases m = 0, 1. In the considered 
example of formula (|2,9j) one gets the profile (Fig. ) for the long arc 7* and the profile 
(Fig. 2) for the short arc jf. 

Let us present the formula for the wave amplitude in a neighborhood of the front 
but outside of some neighborhood of the focal points. As we have just seen in the 
previous example, points of self- intersection can appear for t > t CT . The amplitude of 
the wave in a point x belonging to a neighborhood of these points now is the sum of 
the contributions coming from different ij)j(x, t), yj(x, t), and Sj(x, t) with index j, and 
with the Maslov index m(ipj(x,t),t). 

Theorem 1. In a neighborhood of the front but outside of some neighborhood of the 
focal points the wave field is the sum of the fields 

(2.9) 

Outside this neighborhood of the front 7^ i](x,t) = 0([i 3 / 2 ). Again the function 
F(z,n{ij))) is determined in (|1.17j) . 

Let us emphasize that the number m has a pure topological and geometrical char- 
acter and can be calculated without any relation with the asymptotic formulas for the 
wave field. From the theorem |^1 it follows that, in order to construct the wave field 
at some time t and in a point x, one has to know only the initial values rj\t = o and 
r]t\t=o and has not to know the wave field 77 for all previous time between and t. The 
trajectories and the Maslov (Morse) index take into account all metamorphosis of the 
wave field during the evolution from zero time until time t. 

Remark 2. Finally let us note that it is possible to define the Maslov index in the 
singular (focal) points also (see [?]), actually this definition is associated with a chain of 
covering maps of the front and gives the correct wave field asymptotic in a neighborhood 
of the focal point. It is useful to distinguish these two types of Maslov indices; we also 
meet in our calculations the second one, denoted by m, but after the discussion of the 
index m (see subsection \3.2\) ). 
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2.4 Wave field asymptotic in a neighborhood of focal 
point 

2.4.1 Completely nongenerated focal points and coordinate system 

Now we consider the situation when for some t the point (P F ,X F ) = 
(P(ifj F (t), t), X (i(; F (t) , t)) corresponding to the angle i^ F (t) is a focal one. In this point 
X^=0 and one has to use another asymptotic representation for the solution. Roughly 
speaking the neighborhood of the point X(ip F (t), t) on the plane can include several 
arcs of 7t with the angles ip far from ip F (t). This means that one has to take into ac- 
count contribution of all of these arks in the final formulas for rj in the neighborhood of 
the point x = X(ip F (t), t). The influence of nonsingular points are defined by formula 
()2.9j) and the influence of the points from the neighborhood of the focal points are 
described by formulas (|2,16|) given below. Thus it is necessary to innumerate the focal 
points with the closed projection and write P(ip F (t), i), X(ip F (t), t). These points have 
the closed position X F = X(ip F (t),t), but different momentum P F = P{ipJ (t),t). To 
simplify the notation we discuss here the influence into rj only one focal point omitting 
subindex j but keeping P F . 

We present the corresponding formula under the assumption that some derivative. 



C F = C(X F ), X F = X(^ F (t),t) = P F = P^ F (t),t), 

■J, - del ( -V r , P F ) = ^TT^ , J ( f> = det(X F ,X^ , (2. 12) 



and the derivatives X^ ' = for 1 < k < n. It means that this focal point is not 
completely degenerate. For future it is convenient to introduce the "mixed" Jacobian 

J = det(X,P^,t) = ^^(^) (2.11) 

and some characteristics of the focal point (P F ,X F ): 

P F C \ 

-F zjF\ _ C F dgt(P, Py,) T (n) _ a „ +( yF Y (n)F, 

Co 

Again the topological characteristic appears, i.e. the Maslov index of this focal 
point or its neighborhood (it is the same), but now it depends on the choice of the 
coordinates in the neighborhood of (P F ,X F ). It is natural to choose the new coordi- 
nates (x'i, x' 2 ) associated with the nonzero vector X F = X(ifj F (t), t); namely we assume 
that the direction of the new vertical axis x' 2 coincides with the vector X F . We put 
k 2 = *(fc 21 ,fc 22 ) = X F /\X F \=X F /C F = P F C F /C , kx = *(fc n ,fc 12 ) = (k 22 ,-k 2l ) and 
introduce the new coordinates p',x' in the neighborhood of (P , X F ) in the phase 
space M.p X by formulas: 

x[ = ( kl , x - X F ) = - det(X V" XF) = -^det(P- x - X F ), 
x' 2 = (k 2 , x - X F ) = ( XF ^ ~ XF ) = Cf {pFj x _ x F h 

pi = (ki,p), p' 2 = (k 2 ,p). (2.13) 



Easy to see that 



det d xfj =Jf - (2 - 14) 
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2.4.2 Maslov index of the focal point. 

As the determinant J 7^ in the focal point (P , X F ), hence the same inequality 
takes place in some its neighborhood, thus J has a constant sign. On the contrary the 
Jacobian J changes sign in this neighborhood. We define the Maslov index m(P F , X F ) 
of the non (completely) degenerate focal point (P F ,X F ) = (P, X)(ip F (t), t) as the index 
m(P, X)(tp, t) of a regular point (P, X) = (P, X)(tp, t) in the neighborhood of (P F , X F ) 
such that the signs of the determinants J{ip,t) and J(ifj,t) coincide. For instance one 
can choose ip = i/j F (t),t = t ± 5, where delta is small enough. This means that we 
compare the sign of J with the sign of J on the trajectory (P, X) crossing the curve 
Tt in the focal point (P F ,X F ) before and after this crossing. 

Example 4. Let us illustrate this definition for the focal points of the example with 
the axial symmetric bank (|2.7|) (see Fig...). We saw that before crossing the focal points 
J > 0, and after crossing J < and the index of the points on this trajectory is 1. Let 
us find the sign of Jf in the focal point. To this end we differentiate the integrals of 
motion (|2,8[) of the Hamiltonian system with respect to the angle ij) and put tjj = tp F (t). 
We find 

(P F ,P F )=0, X F P F p -X F P F p =p sin i; F (t) for rf> = if> F (t). 
The solution of this equation is 

and ^det(X F ,P^) = det(P F ,P^) = - { p° P ^ } sin^ F (t). Obviously the angle ip F \t) 
belongs to the interval (0,7r) and (P F ,X F ) > 0, hence Jp < and m(P F , X F ) = 1. 

2.4.3 The model functions and the wave profile in neighborhood of 
the focal point. 

Now we present the formulas for the wave field in the neighborhood of a focal point 
x = X F . Let us put a = sign(Jp Jp ) and introduce the function (or more precisely 
the linear operator acting to the source function V (1/1,1/2)) 

/oo roo £n+l 

d£. \ pdpV(pn(ip)) exp{ip(z 2 - £zi - a- — — — )} = 
-oo Jo (n + l)\ J 

/oo noo _ £Tl+l 

di ^dpf(pn(^))exp{ip(z 2 -£zi-aj^-—^)}. (2.15) 



We put 



x\ J F _ det(P F ,x-X F ) J F F x' 2 C _(P F ,x 



\ j F 4 ] \—^C^ CoCfF 1 ^ \JfJF\- ' » Cf 

>F vF\ 



Theorem 3. Each focal point (P , X ) on the front jt gives in its neighborhood the 
following contribution 



F 



V 



n-l 

^ V Co\Jf\^ Re[e -if m (P^^) g g (z F )Z F^ F)] + Q( )} (216) 

into the asymptotic of solution rj. If the several arcs of belong to the neighborhood of 
the point x, then one need to summarize all corresponding functions (|2.16|) and (|2.9|) . 
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3 The geometric base of asymptotics: La- 
grangian manifolds, the Maslov and Morse in- 
dices. 

The aim of the next section is to prove Theorems 1-3. But let us first recall the 
geometrical objects and the important properties of the Hamiltonian system (|2.4[) . 
giving an uniform asymptotic solution to problem (|2.4[) including the behavior in a 
neighborhood of focal points, initial moment of time, calculation of the Maslov and 
Morse indices etc. The majority of these constructions and properties are well known, 
we present them in the most simple form and collect them in our paper for giving a 
self-contained treatment. An exhaustive description of the wave fronts and the focal 
points, their connection with the ray method and the semiclassical asymptotic, can can 
be found for instance in mentioned above monographs and papers. MaslovFed, Berry. 
There exist different equivalent definitions of the Maslov index; one of the aims of the 
next subsection is to recall the definition from [?.I37|I4"T] which, in our opinion, is more 
suitable for concrete calculations. 

3.1 Lagrangian manifolds ("bands") and their properties. 

As we have already said, taking into account the fact that after the appearance of 
the focal points the front line can intersect itself, it is convenient to add to the point 
x = X(ip,t) the corresponding momentum component p = P(ip,t), and consider the 
point r = r(tp,t) = (P(ip, t), X(ip, t)) in the 4 dimensional phase space Rp ia ,. Each 
point r(ip, t) is completely defined by its coordinates, which are the angle tp (defined 
mod 2ir) and the "propper time" t. 

Fixing the angle ip we obtain the trajectories (bi-characteristics) of Hamiltonian 
system (j2.4|) in the phase space Rp^, and, fixing the time t, we obtain the front Tt 
in the in the phase space Rp ja; . The projections of the trajectories from M£ x to the 
configuration space (plane) R 2 are the rays. The projection of the curve T t from Rp 
to the configuration space (plane) R 2 are the fronts 74. Different points r(tpj,t) on 
Tt can have the same projection x = X(ipj,t) on 74, but now we distinguish them by 
different angles ipj. 

Let us fix some small but independent of \x number 5. According to [?, ?] changing 
both the angle tp and the time r G (t — 8, t + S) on the cylinder § x (t — 6, t + 5) we 
obtain the 2-D Lagrangian manifold (with the boundary) Mf = {p = P(ip,T),x = 
X(tp,T)\ip G §, r G (t — S,t + 5)}; the angle ifi from the unit circle S and the time t 
from (t — 5, t + 5) G R are the coordinates on Mf, sometimes we shall use the notation 
a = t — t instead of the time t. Actually the manifold Mf has a structure of a 
cylindrical "band" (or closed strip) with the width 25, thus we call it Lagrangian band; 
of course it depends on 5, we omit this dependence to simplify the notation. The family 
of Lagrangian bands M t 2 is invariant with respect to the phase flow generated by 
the system (|2.4j) . This means that the point r(ipj,r) from Mf shifted by the action 
of the flow g^ gives again the point r(tpj,T + t) on Mf 0+t but with the shifted time 
r + 1. Due to definition the coordinate a does not change. That is why the coordinate 
r (corresponding to t) is called the proper time. Sometimes it is possible to choose 5 
arbitrary large, even infinity (e.g. in the case C = const). But in many situation the 
set {p = P(ip,r),x = X(ip,T)\ip G S, r G —00} has the intersection points (e.g. if the 
trajectories P(^,t),x = X(if),r) belong to the Liouville tori), and this set is not even 
the manifold. But for our purpose it is enough to work with the "Lagrangian band" 
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Mf only. Along with the general properties of Lagrangian manifolds, the band Mf has 
very useful additional ones. Let us present them all for the completeness. 
Let us introduce the matrices 

. It is easy to see also that each column- vector ( • J , ( ) and ( ^) satisfies the 



variational system 

5x =Hpp6p + Hp X 5x, 5p = -(TC xp 5p + Ti xx 5x) (3.1) 

Easy to check that these vectors are linear independent and obviously two first vectors 
are tangent to M 2 . 



l)the rank of the matrix ( ^ ) is egua/ to 5 which actually means that dimension of 



Lemma 1. (see e.g./?, 1]) The following properties are true: 

l)tht 
Ml is 2. 

2fBC = t CB which means that M 2 is Lagrangian, 

3)for any positive e the matrix C — ieB is not degenerate. 

Proof. The first two propositions follow from the properties of the variational 
system. It is easy to check them for t = because B = (— VC(0), njj, C = (C(O)n, 0) 
where nj_ = *(— sin^, cos ip). In this argument we use the definition of the trajectories 
(P,X), namely the property P\t=o = n (V ! ))-^k=o = 0, n = 4 (cos tp, simp). Thus 

according to the variational system 1)3.1)1 the vector columns ( ^ ) and ( ^ ) are 



K X J 

linearly independent for each t which gives 1). Also a simple calculation based on the 
variational system 1)3.1 j) gives that j-^BC - l CB) = 0, which gives 2). To prove 3) 
assume that C — ieB is degenerate, then there exists a 2-D vector £ ^ such that 
C£ = ieB(,. Consider the (complex) scalar product 0= 

< £, (*tfC - *C£)£ > c =< 0£, C£ > c - < C£, Bi > c = i{e < C£, C£ > c +- £ < B£, B£ > c ) = 0. 

From this equation it follows that both £>£ = 0,C£ = which contradicts 1). □ 
The same consideration allows one to obtain the following closed result. 

Lemma 2. The propositions of the previous Lemma concerning the matrices B,C is 
true if one changes the matrix B by the matrix 

B = (P-XP,P 4 ,), 

where A = (C x (0), n(ip)). 

Let us recall that the points x = X(i[> F ,t) = X F on M 2 where the Jacobian 

J = detC = det(X,X^) = 

are the focal points 1 . Since the manifold M 2 is generated by the curves Ft as well as 
by the trajectories (P, X) each focal point of one of theses objects simultaneously is 



^^Note that using the Hamiltonian system we can change X by P in last formula as well as in many 
formulas containing X. 
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a focal point for the other ones. Little later we shall show that this definition of the 
focal points coincides with the definition, based on equality =, used in the previous 
sections. 

Let us fix some time t and consider the smooth curve Tt = {p = P(tp,t), x = 
X(ip,t)} on G R£ )X (the "time cut" of M 4 2 ). Then obviously the front j t = {x = 
X(ip, t)} is nothing but the projection of Tt to M. 2 . Hence the focal points on the front 
are also the focal points of the manifold M 2 , and from this point of view the caustics 
of M 2 are called space-time ones. 

Lemma 3. The vector- functions X and Aw, as well as vector- functions P and Aw, are 
orthogonal: (A, A^) = (P, X$) = 0. 

Proof. According to system (|2.4|) the vectors P and X are parallel and it is enough 
to prove the second equality. Let us differentiate (P,X^) along the trajectories of the 
system 1)2.4(1 . We have 

d C 2 BC 2 1 

-<P, = (P, Aw,) + <P, A» = using E3I= —\P\(C X ,X^) + ^r(^^) + ^Tq^' 

BC 1 d(C 2 P 2 ) C\P\ dC _ BC 1 B(C 2 ) BC _ 
1 1 dtp 2C Btp C 1 1 dtp 1 1 dtp 2C dtp 1 1 dtp 

But X\t=Q = 0, thus (P, X^)\t=o = and Lemma is proved. □ 

Corollary. 1) The following equality is true J = det(X, Xw,) = ±|X| |Xw,|. 2) The point 
x = X(tp,t) on the front jt, or the point r = (p = P(tp,t),x = X(ip,t) on T t € M 2 is 
a focal one if and only if J = det(X, Xw,) = 0. 

According to the equality \X\ = C(X) J = det(X, Xw,) as well as the Jacobian in 
some neighborhood of jt can be equal to zero if and only if Aw, = 0. Thus the last 
equation really determines the focal points from the point of view of the Lagrangian 
manifold also. 

Lemma 4. In the focal point x = x F = X(tp F ,t) 1) (P F ,P^) = 0, but 2) P F + 0, 3) 
% = g|det(X F ,P F ) ; where as it was before C = C(0) and C F = C(X F ). 

PROOF. According to the conservation law ((2.5(1 (P, P^)(tp F ,t) = 
c 2 p 

(V( Qzfxj),X^)(tp >*) — 0- To prove the second inequality one can mention that 
the vector-function (P 1 p,X 1 p) T satisfies the linear (variational) system with non-zero 
initial condition. Thus both components of the solution cannot be equal to zero. 
To prove 3) we write ^-|^ = ^f = [det(X,Xw,) + det(X,X^)]^ = ^,F = using 12.51 = 

det(X,P < V^,Xw, >)|^ F+ det(X,(^g^))|^ F = ^det(X,P^)|^. □ 

Corollary. In the focal point 1) % = §|det(X,Pw,) = ±|^^X| |P^|(^ F ,t) / 0; 
2) during the passage through the focal point the Jacobian J changes its sign from - to 
+ if det(X, Pwj)L = w,f > and from + to - if det(X, P^)L = ^f < 0; 3) There exists 
t cr such that J(tp, t) > for < t < t cr . 

Proof. To prove 3) it is enough to note that det(X, Pw,)| t =o = 
C(0)det(n(^),n(^)_i_) = C(0). □ 
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3.2 The Maslov and Morse indices. 



As we said before the front jt as well as the curve Tj is partitioned into arcs with 
the focal points at their ends and it is formed by the ends of trajectories having the 
same topological structure. This means that they have similar crossing (on Mf ) 
with the focal points and the same topological characteristic, i.e. the Maslov index. 
But it coincides with the the Morse index for the considered situation (see subsection 
(|2.3.1|0 . Let us prove this proposition. Let us remind some necessary definitions 
and constructions. It is needless to say that there exist different definitions of the 
Maslov index. The original definition [?] is based on calculation of indices of inertia of 

matrices gg^ | M 2, q^ x \'*1] \m?> a^upf) etc ' Jt is not very convement in practical 
calculation. Thus we want to present below one [?, ?, ?] which, from our point of 
view, is more progmatic for computer calculations. We already pointed out that the 
Maslov index of the points on x G "ft is the index of the nonsingular point r(ijj, t) = 
(P(ip,t), X(tp,t) on the Lagrangian band Mf. According to the procedure from [?, ?, 
?, ?] one needs to fix the index m° in some marked nonsingular point p = P(ipo, C)> X = 
X(ipQ,() on Mq and then to find the change of the argument of the determinant of 
the 2x2 matrix Ce 1 ' 2 '' = (C — ieB) = (X — ieP, X^ — ieP^p) along one of the paths 
described below joining the marked point p = P(ipo, C)> x = X(ip , () with the given 
nonsingular point p = P(tfj,t),x = X(ip,t), more precisely 

m(ip,t) = m(ip ,t ) + Am, Am = - lim^ Arg det (X - ieP,X^ - ieP^) (tp, t)|^ to . 

(3.2) 

From definition (|3.2j) it follows the fact that we used before: the index can change 
(jump) only crossing a focal point. In fact, if the point (V^v) ^ s a re g u l ar point then 
det(X, X^j) is different from zero so the increment of the argument of the determinant 
goes to zero when e goes to zero, otherwise, if the determinant of (X,X^p) = 0, as it 
happens in a focal point, then the increment of the determinant in 1)3.2(1 is different 
from zero when e goes to zero. We know (see Corollary from Lemma 4) that the 
Jacobian J = det C(ip,Q > for small enough positive (. Thus all the points on 
the front 7^ are nonsingular. So we choose one of the point (P(^q, £), X(ipo, £)) and 
put Maslov index m(V>cb C) = 0. It is possible (and natural) to use one of the of the 
following paths. 1) To move first along the trajectory P(ipo, t), A"(V>o, t ) starting from 
r = C until t = t, then to move from the point P(tpo, t), X(ipo, t) along the curve T t to 
the point P(ip, t),X(^, t) changing the angle ■0. As we will see a little later this choice 
is not very convenient from the point of view of computer realization. The choice 
2) is to move first from the point with the angle tpo to the point with the angle tp 
along the curve T^, and then to move along the trajectory with the angle 0, changing 
time from £ to t. It could happen that during the motion along some closed path on 
Mf one can get nontrivial increment of argument of the determinant of the matrix 
C^ 1 ' 2 ^ = (C — ieB) = (X — ieP, X^ — ieP^) and nontrivial Maslov index of this closed 
path. The following lemma shows that it is not so. 

Lemma 5. The increment of argument of the determinant of the matrix = 
(C — ieB) = (X — ieP, X^ — ieP^p) along any closed path on Mf is equal to zero. Thus 
the Maslov index of any closed path on Mf is also equal to zero. 

Proof First let us show that Lemma is true for t = 0. Obviously all closed paths 
on Mq are homotopic one to another. Let us choose as a path the curve r^, 5 > C > 0. 
According to Corollary from Lemma 4 the Jacobian J = detC(ip,Q > for £ small 
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enough. Thus the increment of argument of determinant of matrix C — ieC over is 
zero and Maslov index of this path is also zero. Hence according to the property of 
Maslov index it is equal to zero for any other closed path on Mq. Any path on Mf 
could be obtain from some path from Mq by means of the canonical transform (the 
flow) g\j. But this transform preserves the Maslov index of closed path. □ 

Remark 3. Thanks to this statement the Bohr- Sommerf eld quantization rule does not 
appear in asymptotic solutions, although the paths Yt on Mf are not contractible. This 
fact is also natural for Lagrangian manifolds associated with the asymptotics of Green 
functions for evolution equations in Euclidian space (like van Vleck formula for the 
nonstationary Schrddinger equation). 

Taking into account this Lemma, choosing the the second way for the calculation 
of the Maslov index and putting m(ipo, Q = for some small positive C < 8 and some 
fixed fixed ipo, (and hence for any ip) we can write for the index m(ip,t): 

m(ip,t) = Am, Am = - lim Arg det(X - ieP, - ieP^)(^,t)|J'J. (3.3) 

The Corollary of Lemma 3 allows us (in the considered problem) to simplify the defi- 
nition ()3.2|) . Let us analyze the determinant 

det(X - ieP,X^ - ieP^)(t,^) = det(X,X^) -iedet(X,P^) - iedet(P, X^) - e 2 det(P,P^,). 

(3.4) 

Since det(X, P^) is not equal to zero and the term det(P,X^) = in the focal 
point, then the third term of equation (|3,4|) can be omitted. Similarly taking into 
account the fact that the jump of the index is an integer number it easy to show that 
the term e 2 det(P,P^) does not play any role in the calculation of Am and also can 
be omitted. Thus instead of the determinant det(X — ieP, X^, — ieP^) (t, tp) we can use 
the determinant det(X, X^, — ieP^,)(^,t). 

Thus one needs to found the jumps of m during the crossing of the focal points and 
find 

Am = i £ lim o Argdet(X,X^ -ieP^)(^,r)|;:^+^, (3.5) 

where At > is small enough, and tw, is the time at which the trajectory crosses the 
focal point (with coordinates ^,£</,)with the angle ip. There may exist several such t^p, 
but all of them, according to the point 1 of the Corollary of Lemma 3, are isolated 
with respect to t. In fact the derivative of the Jacobian is different from in a focal 
point so the Jacobian has an isolated zero and so those zeros cannot accumulate in one 
point. Now we use again this Corollary. Obviously we can change r with t F in the 
right hand side of (|3.5|) . But the term det(X, P^,) (?/> F , t F ) characterizes the increasing 
or decreasing of the first term. Hence if det(X, P^,) (^ F , t F ) < 0, then the argument of 
the complex vector in (|3.5j) changes on the upper half plane from O(e) to ir + 0(e) and 
Am = 1. If det(X, P^,) (ip F , t F ) > then the argument of the vector in (|3,5|) changes 
from 7r + 0(e) to 2tt + 0(e) and again Am = 1. Thus we obtain the following important 
result. 

Lemma 6. The Maslov index m(tp,t) of any nonsingular point (p = P(ip,t),x = 
X(tf),t)) G Ft with the projection x = X(tp,t) on the front is equal to number of focal 
points on the trajectory P(ip,Q,X(ip,()iC £ *- e - it coincides with the Morse 

index of this trajectory. 
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3.3 The behavior of the front near the focal points. 



For the future developments it is useful to have the description of the wave front 7$ in 
a neighborhood of the focal points. The focal points which are ends of arcs 7? of the 
wavefront belong to the caustics, a well known concept in geometrical optics and in 
the space-time wave theory. It is an important fact that in our problem the caustics 
do not depend on the time because the family of manifolds (bands) Mj* is invariant 
with respect to the phase shift gh. One may distinguish two types of sets organized 
by the focal points with stable and unstable structures with respect to small changes 
of the Lagrangian manifold M f 2 (or of functions X(ip, t), P(ip, t)). The first type is 
under the so-called "general position", there exist only a finite numbers of them and 
in the considered 2-D situation namely there are only two types: the so-called fold and 
cusp[?]. These curves are presented in the Fig.l. Sometimes there exist also different 
focal sets with unstable structure. For instance the circle Tq = p = n(ip),x = on 
Mq, has the point x = as the projection Tq from Mj, (or Rp X ) to M.^.. Rotating a 
little the coordinate system in Mz ^ x one can obtain a small ellipse on instead of the 
point x = 0. But we show below that namely this "unstable singular" circle determine 
localized functions in the asymptotic constructions. 

Taking into account the smoothness of the vector-functions X(ip,t) one can easily 
describe the behavior of the wavefront near the focal points. Namely let us fix the 
time t and let tJ) F = ip F (t) define the angle (coordinate on Tt) of the focal point 
X F = X(ifj F ,t). We put y = Atp = tp — tp F . Let n > 2 be the minimum degree of 
the Taylor expansion of the function X around ip with increment y. We say that the 
focal point X F is not completely degenerate if n 7= 00. From the point of view of this 
definition the point p = Po = n(ip),x = Xq = is a complete degenerate one. 

Shifting the origin into the focal point and rotating the coordinates we can "kill" 
the second component of the n-th derivative vector X^\ip F (t) , t) and write x[ = 
ay n + 0(y n+1 ), x' 2 = by k + 0(y k+1 ). Here o / 0, 6 / are Taylor coefficients, the 
integer k > n and the prime indicates the new coordinates; actually a and b depend 
on the time t, but now for us it is not important. It is clear that the previous Lemmas 
are true also in the new coordinates. 

Lemma 7. In the non degenerate case only one opportunity is possible: k = n + 1. 

Proof. In new coordinates 

, _ (ay n + 0(y n+1 ))\ , _ fany^ 1 + 0(y n )\ 
\by k + 0(y k+1 ) J ' ^ \bky k - 1 + 0{y k ) ) ' 

Since the vector P' is orthogonal to X'^ everywhere, and, according to the conservation 
of the Hamiltonian H, P' 7^ 0, then we can write 

, _ f-bky k ~ n + 0(y k - n+1 ))\ , _ f-bk(k - n)y k - n ~ l + 0(y k ~ n ))\ 

P - q { an + ay + 0(yi) ) and P + ~ Q { a + 0(y) J' 

where the factor q(t) / 0, qa(t) is proportional to the Taylor coefficient just after 
qa{t). Taking into account point 1) of Lemma 4 we immediately find a(t) = 0. But 
from this it follows that if n, k 7^ 00 and k > n + 1 then in the focal point P^ = 
which contradicts point 2) of Lemma 4. □. 

Corollary. Let n / 00 then in the neighborhood of the focal point x = X F in new 
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coordinates 



,_( ay n + 0{y n+1 )) \ , _ ( any"- 1 + 0(y n ) \ 
\b y n+l + 0{y n + 2 ) ) ' A * \b{n + l)y n + 0(y n+1 ) J ' 

_ (-b{n + l)y + 0(y 2 ))\ , _ ( -b(n + 1) + 0(y))\ 

Q { an + 0{y 2 ) )> ^~ q { 0(y) ) 



(3.6) 



and 



det(X,X V) ) = _^C(X F )y n - 1 + 0(y n ), det(X,P^) = |q|b(n + l)C F + 0(y), (3.7) 

where it was before C F = C(X F ). The last two equalities do not depend on the choice 
of the coordinates. 

Thus in agreement with the Lemma 3 the determinant det(X, P^,) does not change 
its sign in a neighborhood of a non degenerate focal point. 

Finally in the non completely degenerate case, in the neighborhood of the focal 
point, we have X[ = ay n + 0(y n+1 ), X' 2 = by n+1 + 0(y n+2 ), P{ = -qb(n + l)y + 
0(y 2 ), P' 2 =qan + 0{y 2 ). 

Omitting the higher corrections we find the equation for the part of the front in 
the neighborhood of the focal point x F : x' x = ay n , x' 2 = by n+l . The sign of ab for odd 
n and the sign of a for even n defines the direction of the passage from higher to lower 
leaves for odd n and from left to right leaves for even n, see Fig. 5. Let us note also 
that in the general case n is equal to 2 or 3 only |HZ1 EEj • 

It is convenient to express the coefficients a, b, q via P, X and their derivatives in 
the focal point ij) F [i). Putting in formulas ([3.11)1 y = we find 

p[ =0 , P' 2 = qan, P^ = -bq(n + 1), for ^ = ip F (t), 

and 

x ,(n) x d n X , p x' V> P^{n-l)\ p 

a = — = nl^F> 6 = - (n-l)!(n + l)Pr g= x[ (n) for ^ = ^ ^ 

The directions of the vectors P and X coincide in the focal points. Thus we see that the 
coordinates with index "prime" could be chosen the coordinates introduced in ()2.13|) . 
This gives: P' 2 = \P F \ = P{ = (k^Pj) = - det(X F , P F )/C F = -J F /C F , 

• (n) _ /U v(n)F\ = _ j„ f /yF v(n)F\ i n „ = _ r(n) 



X'to = (k!,X(") F ) = -det(X F ,Xto F )/C F = -JP/C F and 

4 n) , nJ F jP _ Cb(n-l)! 



n\C F (n + l)!C F C ' H j(») 

F 



(3.8) 



3.4 The jumps of the Maslov index along the front. 



Let us find the jumps Am of the Maslov index during the passage through the focal 
points along the front. We fix the time t > and consider the path to cross non 
degenerate focal points (studied above) starting from the angle vpo — 5 and ending at 
the angle ipo + S. 

Lemma 8. The following equalities are true (see Fig. 3): for odd n Am = for even 
n and sign(JpJ^) = ±1 Am = ±1. 
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Proof. Similarly to the proof of Lemma 4, taking into account the inequality 
det(X,P^) / instead of ()3.2j) we can write: 

Am = i e lim Argdet(X,X^-kP^)(^t)|^;i = (3.9) 



± lim Arg [-qaClX^y 11 - 1 + 0(y n ) - i £ (q 2 C(X F )b(n + 1) + 0(y))]|^ = (3.10) 
i li.n , A,g[-v»- 1 - M " ,' ' ' ))\^_ d . (3.11) 



1 

7T E— + 11 

I lim Ar g [-y^-ie( qHn + 1) ^= S 

7T e^+0 a 

We see now that the complex vector-function —y n ~ l — ie{ qb ^ n ^~ 1 ^ ) lies in one half plane 
for even values of n for each y and in one quadrant for odd values of n. So for n even 
one has Am = — 1 if qab > and Am = 1 if qab < while for n odd Am = 0. To 
finish the proof it is enough to take into account formulas (|3.8|) . □ 

Coming back to the original variables we can make the following conclusion. 

Lemma 9. During the motion along the front 7t 

1 ) the Maslov index does not change if the path does not cross the focal points or 
if the Jacobian J = det (X, X^,) does not change the sign after the passage through the 
focal point; 

2) let the Jacobian J = det (X, X^,) change sign after the passage through the focal 
point then Am = 1 if the signs of J = det (X, X^) and J = det (X, P^,) coincide in the 
end of the path and Am = —1 if the signs of J = det (X, X^) and J = det (X, P^) are 
different. 

Example 5. 

Let us check the correspondence of this conclusion with the index m(x £ 7j?) for 
the example with the axial symmetric bank. Obviously J > on the arc 7*. Since the 
trajectories coming to arc 7 2 meet the focal point only one time and J 7^ in this point, 
then J < on the arc 7I . We have shown in the Example (jlj) that J < in the focal 
point. 

Thus Am = 1 if we move from the arc 7/ to the arc 7^ and Am = — 1 if we move in 
the opposite direction. Finally we find m(x G 7I) = 1, which agrees with our previous 
conclusion. In the critical moment of time when the first focal point just appears the 
Jacobian J(vp,t CT ) does not change the sign after crossing the focal point. Thus the 
index m also does not change and equals to in all regular points of the front 7*" . 
So we see the difference in the jump of the Maslov index depends on the direction of 
crossing the focal points on the manifold Mf . 



3.5 Canonical planes in the phase space, nonsingular and 
singular maps. 

To construct the asymptotic solution of problem (|2,2j) . (|2.3j) in the neighborhood of the 
focal points we need additional construction, related with the fronts, maps covering 
Lagrangian bands M t 2 , indices of these maps etc. Let us describe them also briefly, 
using notation introduced in [IT] . 

The 2-D planes with the focal coordinates x^ 1,2 ) = (x%, X2), x*- 1 ' ^ = (xi,p2), x^ 0,2 -* = 
(pi,X2), x( 0, °) = (pi,P2) in the phase space ^t x are called symplectic canonical planes. 
It is convenient to introduce the multi-indices / = (1,2) corresponds to the canonical 
plane (x 1 ,x 2 ), I = (1,0) to (xi,p 2 ), 1= (0,2) to (pi,x 2 ), I =(0,0) to (pi,p 2 )- 2 We 

2 These multi-indices indicate the replacement of the coordinate Xj corresponding to the entry zero of the 
pair (a!,02 by the momentum p 3 (with the same number j). 
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denote also pt 1 ' 2 ) = (p\,p2), p( 1,0 ) = (pi,— X2), p^' 2 ^ = (—X\,P2), p( > > = (—x\,—X2)- 
We call / the index of singularity. It is convenient to mark the canonical plane by the 
corresponding index / and write TZ 2 . 

According to general property of Lagrangian manifold one can cover M t 2 by the 
maps Qj j with the numbers j such that there exist one-to-one map from fiy 7 to its 
projection to the canonical plane M 2 .. This means the following. Along with the 
matrices B^ 1 ' 2 ^ = B, C^ 1 ' 2 ^ = C it is convenient to introduce the matrices 

B^ 2 \^r) = 

B^HM=[ ~, -,). ir.r)- I" (3.12) 






(dPx 


dPi\ 




dt 


dib 


= ( 


8X 2 


dX 2 




\ dt 


dip J 






dX! \ 




/ dt 


dip 




dP 2 


dP 2 




V dt 


dip ■ j 



/Ml /§Pi dPi\ 

B^{^r) = -C = [_ d % _ d % , C <°-°W)=B= 4 8%) (3-13) 

\ dt dtp J \ dt dip 7 

The matrices C 1 give the Jacobians J 1 (if), r) = det C^. Then in each map Qj* (r) 7^ 
0. The maps with the indices Ij = (1,2) are nonsingular ones, all others are singular 
ones with the focal coordinates x 1 '. Note that for practical application sometimes 
it useful to choose some rotated coordinates (x'i,x' 2 ) and [p'^p^ in some maps - J . 
It is important to remember, that the Jacobians J = j( 1,2 ) and j(°>°) are invariant 
with respect to the rotation, but the Jacobians J^' 2 \ j( 1,0 ) are not. Actually in the 
considered problem the maps with index I = (0, 0) are not needed. 

Lemma 10. For any time t there exists a finite covering of the neighborhood ofTt from 
Lagrangian manifold M 2 by the maps flj 3 with the indices Ij = (1, 2), Ij = (1, 0), Ij = 
(0,2). 

Proof. It is enough to prove that at least one of the Jacobians j(°' 2 ) or J^ 1 ' ) is 
not equal to zero in each focal point. Assume that both j( 0,2 ) = and J*- 1 ' ) = in the 
point ip F ,t F . Since = in the focal point, this means that X\P2^p = X^Px^p = 
and det(A",P^=0. But according to Lemmas 2,3 the vector P^ ^ in the focal point 
and P^p is orthogonal to X, which is nonzero everywhere. This contradiction proves 
this lemma. □ 

Remark 4. 1) It is possible to prove a similar proposition about the whole manifold 
M|, but, since Mj, is noncompact, the number of maps is infinite, and we do not 
need this fact. 2) For practical applications sometimes it useful to choose rotated co- 
ordinates (x'i,x' 2 ) and (p'x,^) * n some maps ■ For instance in the neighborhood 
of non degenerate focal points X(ip F (t),t) in the situation considered in subsection 
()3.3|) it is convenient to choose the axis x' 2 coinciding with the vector X(ip F (t),t) (see 
subsectionfocal ) . 

Example 6. 

1) Construct first the covering of the neighborhood of the "initial" curve Tq = {p = 
n(tpj),x = 0} on the Lagrangian manifold M 2 . The projection of Tq from the phase 
space M.p X on the space M 2 is shrunk to the point x = 0, in this point X^ = for each ip. 
Thus according to our classification the point x = is a completely degenerate point of 
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the "front" jq. Using the definition of P(t, ifi) , X (r, ifi) as solution of the system (|2.4I) 
we find for r = t = 0: 

P = -C' = -VC(x)\ x=0 , X = n(ifi)C(x), = n ± (V), ^ = 0, 

where C' = — VCXx^—o is a vector with components C\q = J^-(O), C20 = Jf^(0)- 
Hence J (0 ' 0) |r = J|r=o = det(X, X<p)\ T=0 = 0. Btrf 

J (1 ' 0) |r = det ° .) = Cocos 2 ^, J^|r = det ( = Co Wifi. 

\ -C 2 cos^y \C smifi J 

(3.14) 

T/ms £/ie neighborhood of the curve Tq can be covered by four singular maps (see Fig....): 
the map Q^'°\ covering the arc with the angle — f — C<V ; <f + C) the map Q^' 2 \ 
covering the arc with the angle f — C < ^ < if + C; ^ e ^3 ? covering the 

arc u>i£/i i/ie ang/e — C < "0 < x - ' -1 ' anc ^ ^ e ma ^ ^1 > covering the arc with 
the angle — Q < ifi < -x + C- ^ere and furthermore in these examples ( is a small 
enough positive number. The choice of these maps is not unique but the final results 
does not depend on the particular choice it is necessary that they are 4. 

2) For t less than t CT one can cover the neighborhood ofTt by two nonsingular maps 

(1 2) (1 2) 

(see Fig....) e.g. Q,\ ' covering the arc with the angle — £ < ifi < it + and fi^ ' 

covering the arc with the angle n — Q < ifi < 2n + Q. Actually these maps disappear 

in the final formula for wave field (|2.6j) because both of them are nonsingular; they are 

used like temporary objects. From the other side changing these maps by their union 

gives the tubular neighborhood of Tt which is not a map. 

3) Let us construct the covering of the front Tt corresponding to the example of the 

axial symmetric case and for time t = t CT , when the first focal point has just appeared. 

Then we have only one focal point on Tt cx , the one with the angle ifi F {t cr ) = ir/2 (which 

follows from the symmetry of the problem). Obviously we can choose covering consisting 

(1 2) 

of two maps: nonsingular Q\ ' covering the arc with the angle n/2 + 2£ < ifi < 
2tt + tt/2 — 2( and singular Q^' 2 ^ covering the arc with the angle tt/2 — ( < ifi < 7r/2 + £. 

4 ) Now we describe the covering of the front Tt but for t > t cv . Denote by ifi F (t) the 
angles of the left and right focal points. We can choose the maps in the following way: 
the nonsingular map , covering the arc with the angle j—C<ijj< ifi F (t) — (,, the 
singular map f^ 0,2 ^ , covering the arc with the angle ifi F (t) — 2( > < ifi < ifi F (t) + 2^, the 
nonsingular map Q^' 2 ^ , covering the arc with the angle ifi F (t) + C, < ifi < ifi+(t) — C 
the singular map f^ 0,2 ^, covering the arc with the angle ifi F (t) — 2( i < ifi < ifi F (t) + 2(^, 
and the nonsingular map O,^ , covering the arc with the angle ifi F (t) + C, < ifi < + 

Remark 5. 1) We mentioned that in the singular maps covering the isolated focal point 
with the angle ifi F (t) it is convenient to rotate the coordinates in such a way that the 
new axis x' 2 has the same direction as vector X(ifi F (t)). In the new coordinate system 
the singular map always has a singular index I = (0,2) and focal (rotated) coordinates 
{p'i,x'^). 2) One needs to use the covering from item 3) of the previous example not 
only for the critical time t cr but also times for near times t. 
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3.6 The Maslov index of a singular map. 

The last object we need is the Maslov index of chains of maps {07 (*)}. To find it one 
following matrices ci , Ce 2 ^ with the elements 

defined in ()3.1|) : 



3 

has to fix some nonsingular point r(ip, f) in the corresponding map Qj and construct 
their one of the following matrices ci , Ce 2 ^ with the elements of matrices B, C 



Cp) = C - ieB = ( C r 11 ~ i£ l n C r 12 ~ i£B K l2 \ 

C (l,0) = 



Cn - ieBu C12 - «£#i2 

(C 2 i - ie£>2i) cos?? + (#21 + ieC 2 i) sinr? (C22 - ieB 2 2) cost? + (£> 2 2 + ieCzi) sin 7? 

C (0,2) ; 

(Cn - isjBn) cos 77 + (£>n + ieC\\) suit? (C 12 - ieBu) cost? + (# 12 + ieCyi) suit? 
C21 — ie£>21 C22 — T££>22 

c (o,o) : 

(Cn - ieBu) cost? + (£>n + ieCji) sinr/ (C i2 - ^£#12) cost? + (0 12 + TeCi 2 ) sinr/ 
(C 2 i - ie#2i) cost? + (#21 + ieC 2 i) sinr? (C22 - ieB 2 2) cos?? + (#22 + ^22) sinr? 



These matrices are not generated for any 77 6 [0,7r/2] and any positive e in the maps 
with index (1,0), (0,2) and (0,0) respectively 3 . Obviously C^ =0 = = C - ieB 

and Cg|^=o = C 1 — ieB 1 and these matrices determine a continuous non degenerate 
transition from the matrix cl 1 ' 2 ^ to the matrix C 1 — ieB^'^. The corresponding deter- 
minant JJ^ 1 ' ^ = detC^ 1 ' ^ or Jg ' 1 ^ = detci°' 2 ^ is not equal to zero. Let m(ip,f) be the 
Maslov index of the point r(tp,f). Then the index m(O^) of the map is 



m(nV) = m(^,f) + - lim ArgJ*^ • ( 3 - 15 ) 



This definition does not depend mod 4 on £/te choice of the point (ip, f) in a given 
map. The calculation of the index m(Oy') could be technically complicated even in 
quite simple situations. But taking into account the fact that we can restrict ourselves 
to maps with I = (1, 0) and / = (0, 2) it is possible to simplify the application of this 
formula. 

Lemma 11. One can always find a nonsingular point r(ip,f) in the map m(Oy') with 
Ij = (1,0) or (0,2) such that the sign Jacobian J(ib,f) coincides with the sign of the 
Jacobian J Ij (^,T) in this map. Then the second term in ()3.15j) is equal to zero and 

m(fiy') = m(ip, f). 

As the sign of the Jacobian J 1 ^ does not depend on a point in the thus one can 
find it in any point, for instance in a focal one. 

Proof. It is obvious that the second term in IJM.15J) can be equal to 0, 1 or —1 only. 
Consider for e the Jacobian JJq. A simple calculation gives Jq = J cost? + J 1 sinr?. In 



3 As we mentioned above the objects with singular index (0, 0) are not needed in the considered problem 
we present C^ ' ^ for completeness. 
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the interval [0, 7r/2] this function has no zero if J and J 1 have the same signs and one 
zero in the opposite situation. Hence including the parameter e gives only the rule of 
bypassing of the zero point on the complex plane, and it is not necessary to use this 
rule in the case when J J > Thus one obtains lim £ ^ + o ArgJ^ =0 = if a point 
r(i/>, f) is chosen in the way prescribed in the Lemma. Prom the other side, according 
to Lemma 3 and its Corollary, the existence of the focal point in any focal map means 
the existence of nonsingular points with positive and negative signs of the Jacobian J. 
□ 

Corollary. The index m(Qj j (t)) of the singular map coincides with the index of 
any nonsingular point m(ijj, t) on the front T t where the Jacobians J and J 1 ^ have the 
same sign. 

Let us illustrate this lemma by the calculations of indices of maps from the examples 
considered above. 

Example 7.1) Consider the maps from the example © ■ Let us take the points on the 
maps Q,? ,j = 1,2,3,4 with coordinates ip,r, t > 0. We put m(ip,T) = 0. According 

to Corollary from Lemma 4 J(ip, T ) > 0. From this and (|3.14f) we find m(Qj j )=0 for 
each j = 1, 2, 3, 4. 

Example 8. 2)Consider the maps from 3), 4) of example ([5]). Before let us come 
back to formulas Q3.11|l . As we mentioned in remark (|3.5|) it is convenient to rotate 
the coordinates. Then in the new coordinates in the focal point X F with coordinates 

X' = P'C\X-)/C ^ an y and P> = ("^ + X >) . 

Hence J'(°' 2 ) = abC 2 (X F )q 2 ^n(n + l)/Co and the sign of J'(°> 2 ) coincides with the sign 
of the product ab. For both examples 3.) 4-) the sign is negative (see Fig. ) thus the 
index of the map VL^' 2 ^ in 3 of example and indices of the maps Q^' 2 \ r^ ' 2 ^ in 3 
of example (jBJ) are — 1 . 

3.7 Germs of Lagrangian manifolds and their properties. 

The geometrical objects described above one can meet in many asymptotical problems 
having fast oscillation solution. We see that the solution of considered problem decays 
quit rapidly outside of some neighborhood of the front jf This gives the opportunity 
to use the ideas of boundary layer expansions [12]- [M] an d the "complex germ theory" 
|29| . Their geometrical realization contains in a change Lagrangian band by its 
linearization (germ) near the curve T^. 

Definition 1. For each fix t we call the linear germ corresponding to the manifold Mj, 
a vector fiber bundle in the phase space with base coinciding with the front Tt = (x = 
X(ijj,t),p = P(ip,t)) and fibers generated by the vectors X,P. 

Denote a G M. the coordinate on the bundle (which is a linear analog to a proper 
time r), then we can define a family of manifolds A 2 as a strip in a neighborhood of 
the front F t in the phase space ]Rp j2 . 

A 2 = {p = P(tp, t, a) = P{ip, t) + P(ip, t)a, x = X(V>, t, a) = X(tp, t) + X(tp, t)a}, 

(3.16) 
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where ip S S 1 = [0, 2ir], \a\ < «o are the coordinates in A^. It is easy to see that 
M 2 can be approximated by Af , and that the parameter a is used for linearizing the 
functions X(if), t + a), P(ip, t + a) and that it defines a shift of time near the front IV 
Taking into account this fact it is easy to prove the following proposition. 

Lemma 12. 1) With an error of the order 0(a 2 ) the manifold A 2 is obtained from 
Aq by means of a shift of time t along the trajectories of the phase flow with the 
Hamiltonian 7i = H\p\C(x). 2) For the matrices ~B(ip,t,a) = > C(ip,t,a) = 

the following equalities are true B = B+0(a), C = C + 0(a); *CB = t BC + 0(a), 
where as before B = ; C = ■ The last equality means that the manifold 
(band) A 2 is (almost) Lagrangian mod 0(a), the statement 1) means that it is (almost) 
invariant mod 0(a). 

Proof. Consider the Hamilton equations x = Ti. p ,p = —TL X . We expand the 
derivatives of the Hamiltonian 7i around the point X(t,^), P(t,ip) 

H p = H p (X(i>, t),P(i{j, t)) + n pp aP + 7i px aX + 0(a 2 ) 
H x = H x (X(iP, t),P(iP, t)) + H xp aP + H xx aX + 0(a 2 ) 

and substitute in the Hamilton equations and we get the result. The second proposition 
is obvious and use the variational system for the evolution of (x,p): 

X + aX - (H p (X(ip, t),P(il>, t)) + H pp aP + H px aX) = 0(a 2 ) 
P + aP + H x (X(if>, t),P(rf>, t)) + H xp aP + H xx aX = 0(a 2 ) 

and we get the result. □ 

As the germ A 2 is some approximation of M 2 almost all previous proposition, 
geometrical definition and construction (like Maslov and Morse index) related to the 
band M 2 is true for the band (germ) A 2 . From the other side obviously one does not 
need any additional objects besides the family of curves (fronts in the phase space) 

to construct both the Lagrangian bands M 2 and their germs A 2 . It also follows 
from formulas (|2.6f) . (j2.9f) . (j2.16f) that leading term of the solution n based also only on 
these objects. Nevertheless the proof of (|2.6j) . (|2.9j) . (|2.16j) needs something more, and 
it seems that sometimes technically instead of the germ A^ it is convenient to consider 
also the other germ (fiberbundle) namely 

A 2 = {p = P(^,t,a) = P(ip,t)+(P(ifj,t)-\(i;)P)a, x = X(^,t,a) = X(if), t)+X(if}, t)a} 

(3.17) 

where (see Lemma 2) A = (fg(0),n(V>)). This germ also implies the matrices 
B(ip,t,a) = gg^j, C(V,i,a) = CftM,a) = gg^y. But now B = B + 0(a), where 

B = (P — \P,P^). After analysis of our previous consideration it is possible to prove 
the following statement. 

Lemma 13. All previous proposition concerning the germ A 2 and matrices B,C are 
true for the germ A 2 and matrices B,C. 

4 Geometrical asymptotic solution and the 
Maslov canonical operator 

. The central mathematical result of our paper is the observation that the asymptotic 
solution of the problem (|2.2|) can be represented as an integral over dp of the canonical 
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Maslov operator with "semiclassical" parameter h = l/p, defined on the appropriate 
family of Lagrangian manifolds A 2 , and acting on the function V (jl,19|) defining the 
initial localized perturbation. In some neighborhood of the front line jt the final 
formula has the form: 



*? = Re (\/lf j ( 00 K^(VpV(pn(^))dp + o( / x), (4.1) 

and T) = o{jj) outside of this neighborhood. The initial data (jl,19|) . the representation 
of the asymptotic solution (|2.6[) . (|2.9[) out of the neighborhood of the focal points as 
well as the future representation of the solution in the neighborhood of the focal points 
is only a realization of (|4.1|) in the corresponding domain of R 2 . As we said before 
the integral over parameter p in (|4. 1|) plays a very important role: it implies the 
decay of the function r] outside a neighborhood of the front and in turn allows one to 
simplify the objects and formulas appearing in the construction of the Maslov canonical 
operator. This simplification is based on the mentioned ideas of the "complex germ 
theory" [271 I29| . but in a simpler "boundary layer" version. As we said before from 
the geometrical point of view it means that we can use germs A 2 or A 2 instead of the 
Lagrangian band Mf in ()4.1|) . In next subsections we shall describe the functions and 
other objects determining the operator K 1 ?^. 

4.1 The functions on Lagrangian bands M^. 

a. The action-function. The Lagrangian property allows one to define on the family 
Mf, function s(^,t,a) satisfying the equation ds = (P, dX)\ 



Lemma 14. The phases (ip, r) on Mf is equal to Cqcx. 

Proof. First let us find s on Mq . Then the coordinate a is the proper time: t = a. 
But in this case we have 



r{ip,T) r(0,r) r(i>,r) 

/ (P,dX) = / (P,X)dt + / < 

'(0,0) J{0,0) J(0,t) 



The second term in the last expression is equal to zero according to Lemma 2. Changing 
X by the right hand side from system (|2.4|) and using the integral of motion (|2.5j) . 
we find (using the proper time r) s(ip, r) = Cqt = Cqo. According to HS1 1^| 
to construct the action on the band M t 2 one has to add Jq Cdt to Co a, where the 
Lagrangian C = (p,TC) — TC- But for the wave equation C = 0, which gives the first 
proposition of Lemma. □ 

The phase s implies the phases associated with projection to singular maps fij: 

sWty, r) = C t - P 2 (^, r)X 2 (^, r), s (0 ' 2) (V, r) = C t - P^, r)X x r). 

The initial source function V(y) implies on Mf the function (more precisely the 
family of function depending on parameter p £ (0, oo)) 

fM) = JpV(pn(il>). (4.2) 

We need also the smooth cut off function eo(a), eo(a) = 1 for |a| < ao, and 
e o( a ) = for \a\ > 2aQ where ao is some small enough positive number. The product 
f(p,tp)eo(a) gives the finite function on the band Mq. We continue it on all family 
Mf assuming that it doesnot depend on time t. 
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b. Functions in the maps In each map 0^ the Jacobians J 1 are not equal to 
zero. This means that one can construct the smooth solutions (ipj (sc ), ^(x 1 )) G 
of the system of equations 

X I ^,r)=x I . (4.3) 
Let us emphasize that it could exists several angles ip^x 1 ), corresponding to one vector 
x . Although the manifold Mf is invariant with respect to the shift g^ it does not mean 
that there is no dynamics on Mf. It only means that turning on the time dependence 
we transform the objects related to Mf in a special way. Namely let us use from the 
beginning the coordinate a instead of the proper time r. We have already mentioned 
that the points r(ip, a) G M^ after the action of the transform gt^ passes to the points 
r(tp, a) G Mf with shifted coordinate r = a +t, but with the same angle tp. Thus the 
equations 1)4. 3 j) are changed by the equations 

X I (il),a + t) = x I ; (4.4) 
The following trivial proposition is very useful. 

Lemma 15. Let ipj(x I ),Tj(x I ) be the solution of the equations (|4.3j) in the map (£) 

and let the point r(ip,a + t) G Mf with coordinates ip,a + t belong to the same map fiy 7 . 
Then the angle component ipj of the solution of the equation 1)4.4)1 does not depend on 
t: ipj = ipj(x j ) and atj = ctj(x I: >,t) is 

ipj (x Ij , t) = ipj (x Ij ) , ctj (x Ij , t) = Tj {x Ij ) - t , (4.5) 

Proof is obvious □. 

Using (|4.4j) and Ij4.5|) we can rewrite action-functions and the Jacobians in the 
coordinates x 1 . The behavior of the functions on Mf is different with respect to the 
shift g^. Namely the function s = a, and the functions /,e are constant , this means 
that for each t in the coordinates ip, a one has the same form. On the contrary the 
functions s 1 and all the Jacobians 4 J 1 take the forms: 

s (1,0) (V>, a, t) = C a - P 2 (ip, a + t)X 2 (ip, a + t), 

s {0 ' 2 \ip,a + t) = C a- Pi^a + tyX^a + t), (4.6) 

J 1 = J T (ip, a + t). (4.7) 

Now in the maps {£lj j (t)} we want to pass from coordinates ip,a to coordinates x Ij . 
This gives us the actions, Jacobians etc. in the coordinates x ■? : 

Sf' 2 \xi,x 2 ) =T j (x 1 ,x 2 ) -C t, 

Sj\xi,p 2 ) = Tj{xi,p 2 ) -p 2 X 2 (ipj(xi,p 2 ),Tj(xi,p 2 )) - C t, 

Sf ,2 \pi,x 2 ) = Tj(pi,X 2 ) - p 1 X 1 (ipj(p 1 ,x 2 ),T j (p 1 , x 2 )) - Cot, 

J I /(x I i) = J I i^j(x I i),Tj(x I i)), 

e* = e(T j (x I i)-t). (4.8) 

Let us emphasize that the complicated notations only reflect the situation: each map 
has it own number j and index of singularity Ij. (See examples below). 

Finally we need to introduce the partition of unity with the maps Qj j (t) covering Tf. 
the set of smooth functions ej(tp,t) associated with the covering {Q^ (t)}: suppej (ip) G 

"j i (t),E 3 -e 3 -W = l- 

4 To simplify notation we do not introduce a new symbol for time-shifted Jacobian. 
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4.2 The time-dependent canonical Maslov operator on 
the invariant manifold Mf. 

Now everything is ready to determine the canonical Maslov operator Kvj 2 , acting 
on the function f(p,ip)e(a), which is constant one on the trajectories of system (|2.4|) 
and depending on the parameter h > 0. It means that this function is the same in 
all points P(if),t + a),X(ip,t + a). Let {Q^ } be a covering of the curve T t . Let us 
divide the set of indices {j} into three parts {j(l, 2), j(l, 0), j(0, 2)} corresponding to 
the maps with indices of singularity (1,2), (1,0) and (0,2) respectively. We put 

V(p,x 1 ,x 2 ,t) =K^(fe) = 

iSf - 2) (0:1,3:2 ,t) 

£ e-f"^^))g lP 5 /(p,^)e J -(V>)|^ J ( gl ^)e(T J -(x 1 ,x 2 ,t) - t)+ 
i6j(i,2) y|Jj ' (X1,X 2 )| 

(1 o) j /•+ 00 exp ' ( ^ (1 '° )(:I1 ' P2 ' f)+:r2P2) 
E e~^ m( ^' )-= / fe -/(p,V)e i (V')|^( a!1>pa )e(r i (gi,p a , <)-f)dp2+ 



i(Sr' 3 ^ (pi,a:2,t)+a:iP:0 



^ e -f m( nf 2 ))^^ / 6XP fe =— f(p,i>)eM\^.^ 1)X2) e(r j ip 1 ,x 2 ,t)-t)dp 1 

(4.9) 

where and are defined in (|4.8|) . Now we can construct the asymptotic solution 
T) to the problem (|2.2|) . We put in the last formula h = p/l. 

Theorem 4. 1) For any T independent of p = l/L, the solution r] to the problem 
(|2.2[) in the interval t £ [0, T] has the form: 

I pC() -iir f 00 

V = Vas + o(p), r/ as = y^Re(e— J ^(p,x 1 ,x 2 ))dp. (4.10) 



T/iis asymptotic, apart from terms of the order O(p), does not depend on the choice of 
the covering {O . 3 }, and functions ej, e. 

2) For each time t € [O, T] £/ie function r\ is localized in a neighborhood of the front: 
the function r/ is equal O(p) outside some neighborhood of the front 74. 

Sketch of half of Proof. Using the results [2Si EHl SU EB] one can show 
that the function Q4.9|) is a leading term of some asymptotic solution modO(h k ) 
to original equation (|2.2jl . where k is an arbitrary big by fixed integer number. We 
introduce the smooth cut off function g(y): g(y) = for y < 1/2 and g(y) = 1 for 
y > 1. Multiplying ty k by gip/p 1 ^ 2 ) and integrating the product by dp we obtain that 
the result is asymptotic solution of (|2.2|) modO(^ 2 ). Then as in [2B] we show that the 
influence of the term /^(l — g(p/ p 1 l 2 ))^ k \t=odp into the solution (|2.2|) is o(p), and 
hence the function rj as from is a leading term of some asymptotic solution of (|2.2|) . 
Now we need to check the conditions (j2.3j) . But it is better to do after simplification 
of function and we shall do it in the next subsection. 

4.3 The germ A 2 of the manifold M t 2 and the simplification 
of the asymptotic. 
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Since the function (|4.2[) decays quite rapidly when the point x goes away from the 
front, it is possible to change the functions S Ij , \Jj j in a neighborhood of the front by 
their Taylor expansions. The nice fact is that one does not change the accuracy 0(/i) in 
formula (|4.2ft using only the zero, first terms and sometimes second terms of the Taylor 
expansions of the phases, and zero terms in the other functions. All these expansions 
are expressed via the vector functions (P(ip,t), X(ip,t)) and matrices B(ip,t),C(ijj,t). 

We need the Taylor expansion in the following form. Let the equations (2/1,2/2) = 
(Yi(V'), Yi(ip)) determine a smooth curve T in some domain in R 2 and (2/1, 2/2) be 
a smooth function in some neighborhood T> of T. Let q(ip) be the smooth family of 
nonzero vectors with components (91,92) transversal to the curve T. This means that 
the vectors q{ip) and Y^fy) are not parallel. The parameter ^ on T and the family 
of vectors q(i/j) define the curvilinear system in some neighborhood of T: each point 
2/ in this neighborhood can be characterized by two values: ip(y) and the length (with 
the sign) z = (y — Y(ip(y)) 1 q(tp(y)))/(q^(y))) 2 of the vector y — Y (tp(y)). To find the 
value of tp(y) one has to solve the equation 

(y - y(V0,9±(V0> = (yi - nWOMVO - (2/2 - ^WOMWO = 0. (4.11) 

Lemma 16. The following expansion is valid: 

*(y) = *(y)+(|V), {y -Y))+±((y-Y), ^(Y)(y-Y))\ Y=Y{m) +0((y-Y^(y))) 3 ) 

(4.12) 

Proof follows from the 1-D Taylor expansion of the function ^(Y+qz) with respect 
to variable z □. 

Now we want to apply this lemma to the phases S- j and Jacobians J- 3 in (|4.8|) and 
(|4.2j) . The variable y are x 1 , the curve T = {x 1 = X^(^,t)}, thus the solution ip will 
depend also on time t. We need the first and second derivatives of S 1 in the points 
x = X^ipjt). From the general theory of Hamilton- Jacobi equation [?, ?, ?] it follows 

8 S 1 B 2 f) P 1 

= P*W>,t), ^ = ^OM) = iKfoMXc'CiM))- 1 , 

where the matrices B^C 1 are defined in ()3.13j) . Now let us choose the vector q as 
following. In the case / = (1,2) q = X(ip,t) = (X^)±, then equation (|4,11|) is 
equation. ()4.4|) . In the case J = (0, 2) q = *(0, 1), then equation (|4.11|) is 

PifcM)=pi. (4.13) 

In the case 1= (1,0) q = *(1,0), then Ea. (l4~TTl) is 

W,*)=P2- (4.14) 

We denote ^(x Ij ,t) the solution of these equations in the map with the number j. 5 
Then after some algebra we obtain the following formulas in the maps with numbers j 

J 1 / (x 1 * , t) = {J 1 / (V, t) + 0(x^ - X 1 ! ty, t))}\^ (a , t) , (4.15) 



5 These solutions are different from the solutions ipj(x Ij ) of equation 1|4.5[) . we use almost the same symbol 
to simplify the notation. 
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here 



Sj ll2) t) = (P&, t),(x- Xty, t)) - t)), C x (Xty, t))) (x - Xty, t)\ 

jf* 2) =det{X,X^,t) (4.16) 

sf> 2) ^, t) = -p^, t)x x {$, t) + p 2 (^, t)(x 2 - x 2 y>, t))+ 



-(x 2 -x 2 ^,t)y 



2 P\X 2 ^ — X 2 Piip 

J^H^t) = detC(°' 2 )(^,t) = {P 1 X 2i) -X 2 P Xi) )^,t) (4.17) 
S^°\ip, t) = -P 2 (ip, t)X 2 (^, t) + P 1 &, t){x x - X^ip, t))+ 



1 (x 1 -X 1 ^t)f PlP ^- P2P ^ 



2 X\P 2 ,p — P 2 X\^, 

J^°\iP,t) = detC^ftM) = {X X P^ - P 2 X^)(^,t),t) (4.18) 

Remark 6. It is important that the last formulas do not depend on the choice of the 
vector q with the same accuracy they are valid. 

Theorem 5. The proposition of Theorem 1 is valid if one changes in the formu- 
las g2J,g2J S 1 / by S 1 /, J 1 /, by J- J , Vi(^) by ^'(a^,*), efatfi ,t) - t) by 
e(\x Ij — X T i (ipj (x Ij ,t)\). 2) In the singular maps in formulas (|4.9[) one can change 
the integration over pj £ (—00,00) by the integration over the angle tj) € Hy*, putting 

Qp. 

pj = Pj(ip,t) and dpj = -Q^-(^,t)dip adjusting the limits in the integral with these 
change. 

The idea of Proof. The proof in regular maps is no more but the Taylor 
expansion of regular components in (|4.9j) . (|4.2|) with respect to distance from 7^. The 
proof in the focal maps based on Taylor expansions but also on estimates of some 
rapidly oscillating integrals. 

4.4 Derivation of formulas from Theorems 1-3 

Now let us apply this Theorem in different cases. 

a. A verification of initial data We shall start from t = 0. It was shown in 
Example (jHJ) that one can cover the neighborhood of the curve To by maps fi^ 3 and 

O^i ■ Taking into account the definition of the functions ej the equalities P(ip, 0) = 

n(VO,X(^,0) = 0, P(V,0) = -0,(0) = -C' ,X(V,0) = Q)n(V0 and m(of 0) ) = 



32 



m(fr. ' 2 ') = we easily find in these maps 

<|(/e°)|t=o = 
T/4K x|(C,(0),n(^)) 



f(p,ip)exp (-i : 2 — )ej(^)e (xi)cos^dV- 

,,y4_ s I cos -01 v 2/iC sm^V 

i- /< , *ff^i cos i/' + ^o sin tb) „ 

; f(p, exp H 1J) 2 ; ) e j W e ° xi cos 

37T/4-? I cos^l 2/iCosm^^ 7 

— — / p, V exp (i ^ y ; ' 2 V y 7/ ) e j V e° (x 2 ) sin V# - 

| sin 2/iCo cos z f/j 

/(p^)exp (» a ^ (0) ' n /f ) )e J -We°(x 2 )8inV>#}. 



'57r/4-s I sin ^ | 2/iCo cos 2 ^ 

Discarding the absolute values under the square roots we see that sines and cosines 
in nominators and denominators are cancelled and only signs + appear before the 
integrals: 

. 7r /4_ f 2/iCo sin ^ 

37r/4+? !(£i£M±h£2£Hii) ,, M ,.x 2 (CJ0),n(ib)) N , . .,, 
»r/4-s 2hC cos z ip ' 

e k / p,^ exp H • 2 ; ejWe(xi)#+ 

37T/4-? 2/iCo sm z V 

e k / P,^ exp i^iAW/ e :c 2 # . 4.19 

Now let us put h = p/p and integrate this expression over p from to oo. 
Lemma 17. The following asymptotic equality is true: 

^Re(e"-/ 4 dp VpK^V(pn(^)e)| t=0 ) = V(x/p) + O(p). 



2tt 

Sketch of Proof. Let us first consider the last integrals without the factors 

exp (- ^ 2 2fe^° s in" ( ^ )) ) e ( 3; i) and ex P { ^hcflS^ ) e ( x 2)- Then taking into account 
the equality ^ • = 1 we obtain the integral 



^ /*00 /*27T ip(x-^ cos fp + :E2 sin ^0 



o jo 



(V(pn(i/;))pdilHLp, (4.20) 



which without the factor , 8 _ is no more but the inverse Fourier transform of the 
function V{p). Thus this integral is V(x/p). Now we need to proof that the factors 
exp (-i ^ggW) ) e ( gl ) and exp (j ^^gW) ) e (g a ) change the function Q 
only by 0{yfp). Consider for instance the expression corresponding to the first in- 
tegral in (|1~lTJ)) 

/ / e ■ F(pn(V (1 - exp {-ip 2 \ xK h 2 ^ >' )) ej ^)pal^dp 
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We divide this integral into two parts with the help of the cut off functions e p (p/\fti) 
and 1 — e p (p/Vh). Here e p (z) = 1 for z < 1/2 and e p (z) = for z > 1. It is simple 
to estimate the first part, which gives O(h). The second part we can integrate by part 
which gives the required estimate. The estimates of other integrals are similar. 

The verification of the second condition in (|2.3j) established by differentiation and 
similar calculations. 

Remark 7. The constructed solution could be decomposed into two parts. The first 
part is based on the used construction but with the different function f. The second 
part is similar, but based on characteristics and family of invariant Lagrangian bands 
associated with the Hamiltonian — \p\C(x). But all calculations related with this part is 
equivalent to complex conjugation of the second one, which allows one not go beyond 
to this negative Hamiltonian. 

b. Asymptotics corresponding to regular points. Now consider the neigh- 
borhood of the regular point. Then each component from the first sum in (|4.9|) gives 
the following component in the solution (|4.2|) 



V^dC^ „ , -iw^mm! 1 - 2 )) f°° i-fr, ipS a ' 2) (x l5 x 2 ,t) 



J M 



y/\J^' 2 \M\ 

ej(V»)e*(|x- Jr(V»,t))|^j(a. llXa , t ). 
The Jacobian | jj 1,2) (rp, t)\ is \X\\X^\ = C(X)\X^\, index m(fi( 1,2) ) = m{^{x,t),t) 

S (1 ' 2 \xi,X2,t) 

and the integral can be presented as a function F{— — ^ — '—^{x, t)). The function 

F(z,?p) decays at least as |z| -3 as \z\ — > oo. Thus the solution is localized in the 
neighborhood of the front 7$ and with the accuracy one can omit the cut 

off functions in this neighborhood. Also without disturbing of this accuracy one can 
change Sj ' (x±,X2,t) by its linear part (P(tp,t), [x — X[ip,t)))\^ = ^j ^ x>t y. the quadratic 
correction to the phase changes the solution by 0(/U 3//2 ) only. This gives the proof of 
Theorems 1 and 2. 

c. Asymptotics corresponding to singular maps. Finally let us consider 
briefly and without rigorous estimates the behavior of the solution in the neighborhood 
of a nondegenerated focal point. To this end we fix the time t and the angle ip F (t) 
determining the focal point P F = P(ip F (i),t),X F = X(ip F (t),t) and rotated coordi- 
nates choosing the direction of the new axis x' 2 coinciding with the vector X{tp F ,t F ). 
It gives us opportunity to make all consideration in the map with / = (0, 2) and use 
formulas (|3.8[) . Thus we need to investigate and to simplify the integral 



fm(nf 2) ) » /^exp 



+00 eYn i(s^ 2) (p'iy 2 i)+<p'i) 

f(p, V')ej(V ; )l^=v>,(pi) e *( a; 2 - x 2)dpi 



</j°'Vi,4)l 



Like in the case for t = we change the integration over dp\ by the integration over 
dip = dy. Taking into account formulas H2.13j) we find: 



-Sf ' 2) (pi , 4 , t) + x' x p' x \ v ' l= p{ = 

p[(x[ - x[) + p^ x ' 2 - x> 2 ) + \{ X ' 2 - x>&, t)f p ] p };~ p ^ = 

-qb(n + l)yx\ + qanx' 2 + qaby n+1 + 0(y n+2 ) + 0(y 2 \x'\) + 0(x' 2 - by n+1 ) 2 
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According to $T£M the Jacobian J^°' 2) (p[, x' 2 ) = J F + 0(y). We put dP[ = -qb(n + 
l)dy, which gives in the integral the factor \qb(n + 1)| and the same limits —00,00. 
It is possible to prove that in the neighborhood of the focal point (depended on p) 
after integration over p the corrections denoted as O(-) gives the small correction to 
the leading term of asymptotic of r]. Taking into account this fact and also the fact of 
fast decaying of r/ outside of the neighborhood of the front 74 we can omit the cut-off 
functions ej, e*. Finally we come to the following formula in the neighborhood of the 
focal point X F 

F _ \qb{n + l)|V/uCo 

~ / x 

\J2ir\J F \ 

(0 2) /*°° /*°° 

R e{e-^/ 4 -f m ( n j j ) / dp dyp V( /9 n(^ F ))x (4.21) 

JO J-oo 

exp ip(-qb{n + Vjyx'-^ + qanx' 2 + qaby n+1 ) ^ 

Now we express the coefficients a,b,q via JW, J F ,Cq,Cf using formulas Q3.8JI and 
instead y introduce the variable £ = — \ n+1 y. After some algebra we obtain 
formula lj!D!5|> . 

Remark 8. About the quadratic terms (corrections) in the phases and regularization 
of asymptotic formulas. One can see that actually the quadratic terms disappear in the 
final formulas for the leading term of asymptotics. Nevertheless they play important 
role in our construction. First, without this term it is impossible to show that the 
asymptotic solution satisfy the original equation with necessary accuracy. Second the 
passage from representation acting in the regular maps to the representation acting 
in the singular maps is based on the partial Fourier transform and stationary phase 
method (see \25i \fflp. It is well known that the application of the second one 

based under assumption of nondegeneracy of second derivatives. Thus it is necessary 
to preserve the quadratic terms in proofs, and it is possible to omit the only in final 
formulas. It is also possible to proof that instead of the germ A^ one can choose the 
germ Af. This gives another asymptotic solution, but with the same leading term. 
More over this choice can be simplify the verification of the initial data because the 
exponential function containing x\ and x\ (|4.19|) disappear (there are no quadratic 
with respect to x\ and x<i corrections in the phases) . Thus this choice can be viewed 
as the other regularization, but all of them are natural in such a sense that they based 
on the more precise asymptotic construction. 
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